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Abstract 

Lepton-flavour violating r-decays are predicted in many extensions of the Standard Model at 
a rate observable at future collider experiments. In this article we focus on the decay T—>nnp,, 
which is a promising channel to observe lepton-flavour violation at the Large Hadron Collider LHC. 
We present analytic expressions for the differential decay width derived from a model-independent 
effective Lagrangian with general four-fermion operators, and estimate the experimental acceptance 
for detecting the decay r— >/^/i/2 at the LHC. Specific emphasis is given to decay angular distributions 
and how they can be used to discriminate new physics models. We provide specific predictions for 
various extensions of the Standard Model, including supersymmetric, little Higgs and technicolour 
models. 



1 Introduction 



In the Standard Model (SM) with massless neutrinos, lepton flavour is conserved. However, the current 
neutrino oscillation data indicate non-degenerate neutrinos with large mixing angles [lj, and this, in 
turn, implies lepton-navour violation (LFV) within the SM extended to include massive neutrinos. 
This flavour violation is large in the neutrino sector (where the Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) mixing matrix has large off-diagonal entries, unlike the Cabbibo-Kobayashi-Maskawa (CKM) 
matrix of the quark sector), but very small in the charged-lepton sector, with branching ratios for 
LFV decays suppressed by factors of bm^jm^ [2], and thus well below current and future experimental 
limits. Hence, any experimental signal of charged-lepton-flavour violation would be a clear indication of 
physics beyond the SM. While current bounds from non-collider experiments strongly constrain fi — * e 
transitions, the limits on t — > fi and r — ► e conversion are much less stringent pp. Moreover, as we 
shall discuss in detail below, many extensions of the SM predict LFV in r decays at a rate accessible at 
future e + e~ and hadron colliders. 

In this article we focus on the decay r — > fifxfi which should provide a clean signature to observe LFV 
at the Large Hadron Collider LHC [3j|j In general, the decay r— >eee would have similar characteristics, 
though in many models (including several of the models we consider), the couplings to electrons are 
suppressed relative to the couplings to muons. Moreover, muons provide a far cleaner signal than 
electrons in a hadron collider environment [3]. At the LHC, r leptons are produced predominantly 
from decays of B and D mesons and W and Z bosons. In the low-luminosity phase, corresponding 
to an integrated luminosity of 10 fb _1 per year, one expects approximately 2 x 10 12 and 2 x 10 8 r 
leptons produced per year from heavy meson and weak boson decays, respectively. The r leptons 
from heavy meson decays result in a much softer muon transverse momentum spectrum and are more 
difficult to trigger and analyze. Therefore, in our Monte Carlo studies we only include r leptons from 
W and Z boson decays. With standard acceptance cuts (\rf^\ < 2.5 and pr> > 3 GeV), and requiring 
either two muons with > 7 GeV or a single muon with > 19 GeV for trigger purposes, we 
find acceptances for various BSM models of approximately 25-30%. Thus, even restricting ourselves 
to r leptons from weak boson decays only, and assuming a branching ratio close to the current upper 
limit BR(r fifip) < 1.9 x 10~ 7 p], we can expect approximately 2 x 10 s x 25% x 1.9 x 10 ' ~ 10 
r — > fifj.fi events within the acceptance range of a typical LHC general purpose detector after one year 
of low-luminosity running. With 30 fb _1 of data, it should be possible to probe branching ratios down 
to a level of BR(r — > fififi) ~ 10 -8 [31 [6] at the LHC. [Recent results from B factories also begin to 
test branching ratios down to the level of 5 x 10 -8 [7]; a prospective super B facility [8j could probe 
branching ratios of O(10 -10 ).] 

Note that there is an even tighter experimental limit on the radiative LFV decay r — > fij, with 
LFV can also be probed at a future e + e~ linear collider [4] and at muon or neutrino factories [5]- 



1 



an upper bound on the branching ratio of BR(t —* ^7) < 6.8xl0~ 8 pp. The radiative decay itself is 
difficult to detect at the LHC [;3j, and any model that allows the decay r — > fififl to proceed through 
an effective Tfij vertex only is already constrained by this bound, and will be further suppressed by 
a factor of a from the /2/Z7 vertex. However, models where the LFV r decay can be mediated by the 
exchange of new heavy particles escape the tight bound on radiative transitions. With this in mind, we 
will only investigate models of this type, though we will briefly comment on other models afterwards. 

There are many observables associated with the decay r — > [i\x]±. Since we are focusing on models 
where the virtual particles have large masses M 3> m T , the propagators are replaced by —i/M 2 , and 
the obvious difference between models is the chiral structure of the effective vertices. We choose the 
angle G to be the angle between the polarization of the r lepton and the momentum of the anti-muon 
(assuming t~— >n~ fi~ fi + ) and differentiate models according to their partial decay width in cos 0, having 
integrated over all other kinematic variables. Approximately 85% of the r leptons from electroweak 
gauge bosons will be produced in W decays, which result in characteristic spin patterns. For example, 
the decay W~ — > t~v t produces r leptons with left-handed helicity only, so that the polarization vector 
is antiparallel to the momentum vector and can thus be determined experimentally. We note that the 
polarization of t leptons from D and B mesons or from Z decays cannot be determined in such simple 
manner. The r leptons from these sources should thus not be included in our polarization analysis, but 
instead dilute the signal. 

The paper is organized as follows: first, we will present analytic expressions for the matrix element 
and differential decay width derived from an effective Lagrangian with the most general four-fermion 
(/2//)(/2r) interactions that do not involve derivatives. The result of this calculation has been used in 
a Monte Carlo study to estimate the detector acceptance for various new physics models. In order 
to show how the decay angular distribution can be employed to discriminate between different models 
we provide specific predictions for various supersymmetric, little Higgs and technicolour models, and 
models with doubly charged Higgs bosons. In each case, we will briefly describe the relevant features of 
the model, extract the values for the effective parameters of the general matrix element and differential 
decay width expressions, present the dependence of its partial decay width with respect to cos and 
the acceptance for the decay simulated in a typical LHC general purpose detector with the experimental 
cuts described above. Finally, we will briefly discuss three additional classes of models: those that fit 
the criterion of heavy mediating particles, but have no concrete predictions for the dependence on cos 0, 
those which are very constrained by the r — > ^7 data, and unparticle models. 

2 Model— independent analysis 

We consider an interaction using effective four-fermion vertices. Derivatives have not been included, 
since by the equations of motion these derivatives will be of the order of the lepton masses, which should 



2 



be small in comparison to the scale of the new physics. For other model-independent investigations 
based on effective lepton-flavour-violating interactions, see Refs. [9]. 

The effective Lagrangian is given by (using p and v as Lorentz indices to avoid confusion with the 
symbol /i being used for muon spinors) 

C = G (^g s LL {flP R n)(flP L T) + g s LR {flP R fi){flP R T) + g s RL {flP Lf i)(jlP L T) + g s RR {jlP L p){fiP R T) 

+g V L L {niuPR^l U PLT) + 9 Y R (fl lv P Rf x)(m U PRT) 

+9RL(fiJiy p LV){fif 'Plt) + g\ R {^ u PLp){m v Prt) 

o mJ _ \ / at" _ \ -r / ^ „ \ ( of ^ ^ 



= (prVp a7 V) (pr'PjT) , (i) 

a,6,c 

where -Pt,/,r are the left- and right-handed projection operators, respectively, with Pl/ r = (1 =F 7 5 )/2, 
and a pv = i[y p ,j' / ]/2. The symbols g c ah denote the couplings for the various chiral structures, with c = 
S for scalar, V for vector, T for tensor, so that V s = 1,T = 7^,r T = a pu /*/2, and where a = {L,R} 
is the chirality of the anti-muon and b = {L, R} is the chirality of the r lepton. [The presence of the 7 
matrices on each side of P a in the final line is to allow a to be the chirality of the outgoing anti-muon.] 
The overall constant G is dimensionful, with units of GeV -2 . We will use G to absorb any normalization 
of the couplings g c ah for any particular model, so that we can present the g c ah as integers where possible. 
Note that not all ten couplings are independent, as through Fierz identities one of the four-fermion 
terms can be written as a sum of one or more of the others. However, we choose to keep the extra term 
for convenience. The use of Fierz identities can also show that the omitted g\ L {j l ~^Pr^ {j L ^/2^ >Ll ^) 
and g RR {j^-^^Pl^ (^^TJ-^R 7 ") terms are identically zero. 

Noting that there is a pair of identical fermions in the final state, the transition matrix element is 
given by 

M = GY,9ab{HP^ C l Pal v(PfiHP^ C PMPr) 



a,b,c 



-«(p A4fl )r c 7°i'«7 «(Pp)ti(PMA) rC ^"(Pr)) , (2) 



where A and B label the muons. While we sum over final-state spins when squaring the matrix element, 
we keep the information about the t polarization by using u(p T )u(p T ) = ($ T + m T ) — V where n u 
is the polarization vector of the r lepton. 

Squaring this matrix element and performing integrations over the phase space except for the anti- 
muon energy and the angle between the polarization of the r lepton and the momentum of the 
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anti-muon produces the result for the normalized double-differential decay width: 

1 dr 6x 2 (a + bx + c cos + d x cos 0) 



(3) 



r dx dcos (4a + 36) 

with the approximation that the muons are massless. Here x = 2Ep i /m T is the reduced energy of the 
anti-muon, and the coefficients a, b, c and d are given by 

« = 3| 5 f L | 2 + l2\g v LL \ 2 + 3|^| 2 + 48\g v LR \ 2 + 108|^| 2 + 3\g s RL \ 2 + 48|^ L | 2 + 108|a£ L | 2 + 3\g S RR \ 2 
+l2\g v RR \ 2 - 12 Re {g s LL g v L l + g v RR g s RR ) - 36 Re {g s LR g T LR + g s RL g T R l) , 

b = -2\g s LL \ 2 - 8\gl L \ 2 - 3\g s LR \ 2 - 48\gl R \ 2 ~ 108|oL| 2 " 3|<?£ L | 2 - 48\g% L \ 2 - 108|^ L | 2 - 2\g s RR \ 2 
-8| g v RR \ 2 + 8 Re {g s LL g v L l + g v RR g s RR ) + 36 Re {g s LR g T LR + g s RL g T R l) , 

c = I 9ll 1 2 + 4|5ll! 2 " Mr? ~ 48|<?L| 2 - W8\gl R \ 2 + 3\g s RL \ 2 + 48|^ L | 2 + 108|^ L | 2 - \g s RR \ 2 

-MgU 2 + 36 Re {g s LR g T LR - g s RL g T R l) - 4Re {g s LL g v L l - g v RR g s RR ) , 
d = -2\g s LL \ 2 - 8\gl L \ 2 + 3\g s LR \ 2 + 48\gY R \ 2 + W8\gl R \ 2 - 3\g s RL \ 2 - A8\g^ L \ 2 - 108|^ L | 2 + 2|o^| 2 

+8|o^| 2 - 36 Re {g s LR g T LR - g s RL g T R l) + 8Re {g s LL g v L l - g v RR g s RR ) . (4) 

Hence the normalized differential decay width in cos is 

1 dr If Ac + 3d \ 

- 1 + ^7 cos > ( 5 ) 



T dcos 2 V 4a + 36 

where we use T to denote only the width of the r — > fififi decay, to save on a proliferation of subscripts, 
rather than using it to denote the full r decay width. This normalized differential decay width is now 
independent of the absolute magnitudes of the four-fermion couplings. Muon mass effects are included 
in our full calculation. They are suppressed by powers of m^/rrij- and are thus numerically small. 

3 Model Discrimination 

Here we discuss various new physics models and consider how they may be discriminated using the 
distribution in cos 0. Specifically, we consider the i?-parity-conserving MSSM including see-saw neu- 
trino masses at large tan/3 \10\ lllj. the i?-parity-violating MSSM [12], the Littlest Higgs model with 
T-parity [15] . the topcolour-assisted technicolour model [2], the Higgs triplet model [15], and the 
Zee-Babu model [16 . The decays of r leptons to within each model have been discussed in the 
literature, but, as far as we are aware, only to the extent of predicting total branching ratios in terms 
of the model parameters, except for Ref. .[?] . which considers forward-backward asymmetries for the 
Higgs triplet and Zee-Babu models. 
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(plus diagrams with the muons crossed.) 



Figure 1: Feynman diagrams leading to r — ► fifxfi in the see-saw MSSM with large tan/3. 
3.1 Supersymmetric model with see— saw mechanism 

Neutrino masses can be accommodated in the Minimal Supersymmetric Standard Model (MSSM) by 
adding singlet chiral superfields that have right-handed neutrino components and a large Majorana 
mass term, leading to the supersymmetric version of the "see-saw" mechanism [18\ 119] . However, the 
see-saw MSSM allows for much stronger influence of the mixing in the neutrino sector on the charged 
lepton sector, through large renormalization-group effects in the slepton sector [201 121j . Here, we focus 
on the potential for a large r — > fififl branching ratio escaping the r — > \i r ) bound through an enhanced 
coupling of the MSSM Higgs bosons to muons for large tan /?, as proposed in Ref. [11] 

We start directly from equation (15) in Ref. [11], which gives an effective coupling of the MSSM 
Higgs bosons to a r lepton and a muon. We reproduce it here for ease of reference: 



where in the second line we have explicitly taken the large {a— >/3 — tt/2) and large tan (3 limits. The 
flavour-diagonal muon coupling to MSSM Higgs bosons [22] in these limits is given by 




(6) 



e 



rn 



2-/2 [h° + tan /3H° - itan (3j 5 A°] fx 



(7) 
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Figure 2: 1/T dr/dcos0 for the supersymmetric see-saw model with large tan/3 defined in Ref. 
We display the analytic result obtained for massless muons together with the numerical calculation 
including all mass effects. 



so integrating out the H° and A fields leads to the relevant part of the effective Lagrangian given by 



C cS = -{2G 2 F ) ll4 m T K% 2 {t&nj3f 
= -{2G 2 F ) 1/4 m T K* 32 (tanPY 



m 



sin 9w 2mw 



m 



-tan/? — tan (3j 5 — =- 



m 



ft 



m 



■ —(p-Pl^Prt) 
sm 0w 2rn\v 



(8) 



In the second line of (JSJ) we have used the fact that in the large tjia limit, the mass of the heavier 
CP-even Higgs boson, H°, is approximately equal to m^. We can easily see that the only non-zero 
effective coupling g c ah in the Lagrangian (TT]) is gfjR, which we set to 1. 

Performing the phase-space integrations analytically with massless muons gives the result 



1 dr 



3 + cos 



TdcosG 6 ' ^ 

which can be derived from Eq. and ©. In order to quantify the effect of muon masses, we display 
Eq.([9|) in Fig. [21 along with the result of a numerical calculation incorporating massive muons. As 
expected, the muon mass effects are very small. A Monte-Carlo simulation of the decay r — ► /i/U/i 
within the MSSM see-saw model for a typical LHC general purpose detector (with cuts as described in 
Section [1]) shows an acceptance of 30%. 
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(plus diagrams with the muons crossed.) 
Figure 3: Feynman diagrams leading to r — > in the MSSM without i?-parity. 



3.2 Supersymmetric model with i?— parity violation 

The MSSM without i?-parity (see e.g. Refs. [23] for reviews) includes an additional set of Yukawa 
interactions, between the charged leptons and sneutrinos, which are not necessarily diagonalized when 
the Higgs-lepton Yukawa interactions are diagonalized, allowing for the possibility of charged-lepton- 
flavour-violation by a scalar particle at tree-level. The i?-parity-violating part of the superpotential is 
given by 

= e ab (-XijkLflJjEk + X' ijk LiQ b jD k + KiL°;H^J + -e rs t\" jk U- T))D\ . (10) 

Here, k = 1,2,3 are generation indices, a,b = 1,2 are SU(2) and r,s,t = 1,2,3 are 577(3) indices. 
L, E denote the lepton doublet and singlet left-chiral superfields; Q, U, D denote the quark doublet and 
singlet superfields, respectively. A, A', A" are dimensionless coupling constants and k is a mass mixing 
parameter. 

The terms relevant to the decay r — ► ///i/i are e a b\\ijkL1 L b -E k plus its Hermitian conjugate, which 
give interaction terms ^zHlP-PrT, \* 32 v*lP>Plt , A^^ilM-PrM and X^^il^Pl^- These lead to the 
t — ► /i/U/2 diagrams shown in Fig. [3j Integrating out the sneutrinos, which are constrained to be heavy 
{rriy ^> m T ), gives the following terms in the effective Lagrangian: 



cff 



1 



ID 



-Ai32A*22(P- P KA')(M- P LT) H o~ ^2^*23 {^Pl^^Prt) 



m 



(11) 



from which we can read off that gf L is proportional to Aj32A* 22 an d g^ is proportional to Aj22A| 2 3) 
9lr = Srl = 0' an( l an the vector and tensor couplings are also zero. 

Bottom-up approaches to constraining the MSSM without 72-parity generally consider the minimal 
number of couplings to be non-zero for any individual constraining process [23], e.g. constraints from 
the current upper bound onr-t fifip, assume either Aj32A* 22 ^ and Aj22A* 23 = or the other way 
around, to obtain a conservative upper bound in the absence of destructive interference [25j. Here we 
take two benchmark scenarios: one where gf L = 1 and g^ R = 0, designated "L", and the other where 



7 



Figure 4: Example Feynman diagrams leading to r — > ///i/x in the Littlest Higgs model with T— parity. 
For the full set, see Ref. [27]. 



gf L = and g RR = 1, designated "R". If the decay r — ► fifip, is measurable at the LHC because this 
model is correct, then it is almost certain that it will have been established from other signals, and 
the measurement distinguishing between the two benchmarks would be most useful for measuring the 
i?-parity-violating couplings. 

Performing the phase-space integrations analytically with massless muons gives the result 

( 3-cosG 



i dr 



r dcos © 



for L:g* L = l,g RR = Q, 

(12) 

3 + cos 6 - „ o „ s i 
for R:gt L = 0,g RR = l. 
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Simulation of these decays for a typical LHC general purpose detector shows an acceptance of 25% and 
30% for "L" and "R", respectively. 

3.3 Littlest Higgs model with T— parity 

The additional gauge group(s) and fermion multiplets in the Littlest Higgs model with T-parity (LHT 
model) (see e.g. Refs. [26J) allow for flavour-changing interactions through loops of T-odd particles with 
potentially very different PMNS- or CKM-like matrices. We focus on the way this allows for r — ► n/ifi 
decays as discussed in Ref. [27] . 

We start directly from equations (5.2), (5.3) and (5.4) in Ref. [27J, which give four-point amplitudes 
for T—>fj,fj,fi in the LHT. We reproduce them here for ease of reference: 

Gf e 2 1 -/ 

Al ' = ^^g^ D oddi^)( m r ia ^q P (l+l5))r{p)}[fi(p2h^ (13) 



4 ^|^2^odd D*(Pl)7a(l - 76)T(p)] HP2h°KP3)] ~ (Pl<->P2) 



(14) 



Ahox = ^Cdd[MPi)7a(l-75)r(p)][MP2)7 Q (l-75MP3)] • (15) 

The structure of the amplitude Ay corresponds to setting the couplings g\ L = g^ L = 1 (and all others 
to zero) in the general effective Lagrangian ([1]), while the structure of the amplitude ^4box corresponds 
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to setting g^ L = 1 and all other couplings to zero. The amplitude »4y is already tightly constrained by 
r — > fX'y, so we ignore it in the following. [We also found for input mirror-lepton to mirror-W boson 
mass ratios in the range from 0.1 to 10 that Ay is much smaller than the other amplitudes.] The 
functions Z'JJ^ and Yj^^, defined in Ref. [27], encode all the information from the T-odd particles in 
the decay, and can vary strongly in their relative magnitude. However, we found that for mirror-lepton 
to mirror-W^ mass ratios below 2, Z'JJ^ <C Y'^dd wriUe f° r ratios above about 4 ZjJ^ becomes dominant. 
As a pair of benchmarks, we thus take Z'JJ^ 3> Yj^, so g\ L = g RL = 1 (and all other couplings set to 
zero), designated "Z", and ^odd^^Todd> so g^ L = 1 and all others zero, designated "Y". 
Performing the phase-space integrations analytically with massless muons gives the result 



Simulation of these decays in the LHC environment show an acceptance of 27% for both the scenarios 
"Z" and "Y". 

3.4 Topcolour— assisted technicolour model 

Top colour-assisted technicolour models [Hj , are based on the product of gauge groups SU (3)i x SU (3)2 X 
U(1)yixU(1) Y 2xSU(2) l which is broken to S77(3) QC d x C/(1)em- The group SU(3) 2 xU(l) Y 2 couples 
preferentially to the lighter two generations, while SU(3)ixU(1)yi couples preferentially to the heav- 
iest generation. The generation-dependence of the couplings lead in general to flavour-non-diagonal 
couplings of the heavy Z' vector gauge boson associated with the broken U (1) gauge group once the 
fermions have been diagonalized to the mass-eigenstate basis. 

The top-pion, 7Tj, a condensate which breaks the SU{2)l vacuum, can also mediate flavour-changing 
interactions. Since both the Z' and the top-pion masses and couplings are unknowns, we consider the 
two extreme cases where either of the two particles dominates the decay: we designate the Z' case by "Z" 
and the top-pion case by "P", which occur for example when m nt 3> mz> or m 7Tt <C rnz 1 , respectively. 

We consider r — > ^fifi through the heavy Z' , using couplings as given in Ref. [28]. The flavour- 
diagonal couplings read 




for Z: gY, 



for Y: aY, 





(16) 



FD 

Z' 



-SiCOtfl' Z'^TL^TL + 2f R ^T R ) 

-git&nO'Z'^flL^VL + 2/2 i j7 / >fl + eLl^e L + 2e R ^e R ) 



(17) 



while the flavour-changing couplings are given by 



CW = -\9iZ\[k T ^f L ^^L + 2f Rl ^ R ) + k Te {f L ^e L + 2f Rl ^e R ) 

+ k„ e teLn 2 e'Z'„(fi L ^fx L + 2ji Rl »iiR + e L l^e L + 2e R ^e R )] + h.c. , 



(18) 
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(plus diagrams with the muons crossed.) 
Figure 5: Feynman diagrams leading to r — > in the topcolour-assisted technicolour model. 



where g\ = y/^Ka/cosOw and 9' is the mixing anglj 2 ] for the heavy Z' . The fejijj are flavour mixing 
factors. 

Integrating out the Z' gives the following relevant terms in the effective Lagrangian: 

1 7TQ 



-cff 



MnO'kr,* {fi-fP L HpfP L T + 2jbfP R IJLjjrfP L T 



m? Z t (cos #iy) 2 

+2^P L /i/2 7 ' t P fi T + AJx^Pr^Prt) , (19) 

from which we can read off that g RL = 1, = g^ R = 2, g]^ = 4 and all other couplings are zero. 

Now we consider r — ► //^/i through the 7Tj, using couplings as given in Ref. [29]. The couplings to 
the top-pion are given by the following terms in the Lagrangian: 



-FD 



m t 



+ 



V2F t 
mi 



mt 
+ h.c. , 



(20) 



where m' b ~ 0.1 emj is the part of the bottom quark mass generated by extended technicolour (e <C 1 
is a small parameter). The symbol K denotes flavour mixing factors, v w = v/y/2 where v ~ 246 GeV is 
the scale of electroweak symmetry breaking and Ft is the top-pion decay constant. 

Integrating out the lit gives the following relevant terms in the effective Lagrangian: 



-CcfT 



mimi 

■ -^K* Ti lj{P R - P L )l 3 k(PR ~ Pl)t , 



(21) 



from which we can read off that g^ L = g RR = 1, g RL = g^ R = — 1, and all others are zero. 



S 



2 We note that, in Ref. |28| . the formula for the flavour-changing terms actually has tan 2 # rather than tan 2 #', but we 
assume that this was a typographical error. We also assume that the Hermitian conjugate was also meant to be present in 
the equation. However, we disagree with equation (4) in Ref. [21], where we believe that the numerical factor should be 
7/4096 rather than 25/384 for the case of identical particles in the final state (r— >fifJ,fi, /i/ie, eefi, eee) and 25/1536 for the 
case of no identical particles in the final state (t— tfiep,, flee). 
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Figure 6: Feynman diagram leading to r — > [i\x\i in models with doubly-charged Higgs bosons. 



Performing the phase-space integrations analytically with massless muons gives the result 

7 — 5cos G 



1 dr 



T dcos 9 
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for Z:gy L = l,gl L = 9 y R = 2,gy R = A, 
for P: g s LL = g s RR = 1, g s RL = g s LR = -1 . 



(22) 



Simulation of r — * fj,[i[i decays within the topcolour-assisted technicolour model show an acceptance 
of 29% and 28% for the "Z" and "P" scenarios, respectively. 

3.5 Models with doubly— charged Higgs bosons 

Models with doubly-charged Higgs bosons can mediate LFV r decays through Feynman diagrams like 
that depicted in Fig. El To be specific, we discuss two concrete models, the Higgs Triplet model [15] . and 
the Zee-Babu model [16], both of which have been discussed in the context of LFV r decays previously 
in Ref. [T7]. 

In the Higgs-triplet model neutrinos are given Majorana masses through the addition of a triplet 
of SU(2) with hypercharge 2 (no right-handed neutrinos are introduced) which obtains a vacuum 
expectation value. The doubly-charged components of this triplet can also mediate the decay r — ► /x/i/L 

The Zee-Babu model [16] is another model in which neutrinos are given Majorana masses without 
the introduction of right-handed neutrinos, through the introduction of additional singly- and doubly- 
charged scalars which are singlets of SU{2). The neutrino masses are generated radiatively at the 
two-loop level. 

We start our discussion from equation (61) in Ref. |T7| combined with Table 2 therein. We reproduce 
the (Hermitian conjugate of the) equation and the relevant information from the table here for ease of 
reference: 



cff 



-4G F 

7T 



■}■ 



(23) 
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where 



-4G F n __ J -4G F h^Kp 



_ 53 = and —^54 = , rM (24) 



*2 



for the Higgs triplet model, and 



y " M^ ±± ' ^2 
3 constants i 

charged scalar boson mediating the decay. 



53 = „ , 2 and — 7^54 = (25) 



for the Zee-Babu model. The constants h are mixing factors and Mjr±± is the mass of the doubly- 

L/R 



The structure of Eq. (j23j) to (|25|) corresponds to setting = 1 and all other couplings to zero 
for the Higgs triplet model, while the g^ R = 1 is the only non-zero coupling in the Zee-Babu model. 
We note that LFV r decays with a chirality structure analogous to that of the Higgs triplet model are 
predicted in seesaw models with scalar triplets, see e.g. Ref. [30J. 

Performing the phase-space integrations analytically with massless muons gives the result 

' x _|_ CO s 

for the Higgs triplet model: g RL = 1 , 

2 ' (26) 



i dr 



r dcos G 



for the Zee-Babu model: g\ R = 1 



2 

Simulation of this decay for a typical LHC general purpose detector shows an acceptance of 25% and 
30% for the Higgs triplet and Zee-Babu models, respectively. 

3.6 Other Models 

Finally, let us very briefly discuss some other new physics models, namely those that fit the criterion of 
heavy mediating particles, but have no concrete predictions for the dependence on cosG, those which 
are very constrained by the r — ► /ry data, and unparticle models. 

3.6.1 Heavy particle mediated 

Left-Right Symmetric models [H] introduce a second SU(2) gauge group for the right-handed fermions, 
arranging them in appropriate doublets. The scalar sector is more complex than that of the SM, con- 
sisting of a Higgs bi-doublet and a Higgs triplet for each gauge group. The doubly-charged components 
of the triplets can mediate lepton flavour violation [T7j, but we find that the model is not predictive 
enough to meaningfully differentiate it through 0: the slope of 1/r dr/dcos0 depends on which 
dominates the decay and on the mixing matrices, which are entirely unknown in the right-handed sector. 

3.6.2 Photon mediated 

Any model that predicts lepton flavour violation through loops will predict r — > fi/xp, through a virtual 
photon. However, unless the model predicts other processes that can facilitate r — ► ju/x/i (such as the 
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see-saw MSSM), the bounds from r — > ^7 mean that r — * \jliiix from such models can barely be observed 
at LHC, let alone with sufficient statistics to discriminate models through angular distributions. 

3.6.3 Unparticle mediated 

The phenomenology of "unparticles" [32] has recently attracted a considerable amount of attention. 
In particular, Refs. [33] and [34] consider the effects of unparticles mediating lepton flavour violation. 
However, the hypothesized unparticle sector has no predicted form for the couplings to the SM particles, 
and so has no unique effect on dr/dcos 0. The unparticle nature of r — > fj,fj,p, would instead show up in 
the differential cross-section with respect to the energies of the muons, analogously to the analysis of 
the electron energies in Ref. [34J. 

3.7 Discrimination potential at the LHC 

We summarize our results for the differential decay distributions 1/r dT/ dcos within the various new 
physics models of Section [3] in Fig. [8] and [9] and Tab. HJ [While Fig. [8] and [9] display the numerical 
results including all mass effects, Tab. [T] quotes the analytic results for massless muons, which are an 
excellent approximation to the numerical calculation with full mass dependence.] We can see that 
determining the slope to within 10% should be sufficient to distinguish between the see-saw MSSM, 
the Littlest Higgs model with T-parity and the topcolour-assisted technicolour model. The i?-parity- 
violating MSSM is not distinguishable from the see-saw MSSM for A^AJ^ much larger than Aj32A* 22 , 
but is easily discriminated from the Littlest Higgs model with T-parity and the topcolour-assisted 
technicolour model, since the magnitude of the slope is bounded by 1/6, which is less than the slopes 
of the other models. 

With an estimated 2 x10 s r leptons produced from W bosons in the first year of low-luminosity 
operation of the LHC, the current upper bound on r — > fifj.fi means that even if the branching ratio is 
just below the bound and the full detector efficiency is not much lower than the simulated acceptance, 
we only expect about ten events in the first twelve months. However, after a year of high-luminosity 
operation, we expect as much as a hundred events, which should be sufficient to measure the cos 
dependence of the decay width to within the desired 10%. 

4 Conclusion 

We have analyzed the lepton-flavour violating decay r — ► fififi and presented analytic expressions for 
the differential decay width derived from an effective Lagrangian with general four-fermion interactions. 
The results have been used in a Monte Carlo study to estimate the experimental acceptance of r — > /i/ip, 
decays at the LHC for generic sets of models. We have derived specific predictions for five classes of 
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COS© 



0.5 



Figure 7: 1/T dr/dcos0 for the supersymmetric models with see-saw mechanism (upper figure) and 
-R-parity violation (lower figure) discussed in Section The results are obtained from a numerical 
calculation including all mass effects. 
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Figure 8: 1/T dr/dcosG for the Littlest Higgs model (upper figure) and the top colour-assisted techni- 
colour model (lower figure) discussed in Section [3j The results are obtained from a numerical calculation 
including all mass effects. 
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Figure 9: 1/r dr/dcos© for the models with doubly-charged Higgs bosons discussed in Section [3j The 
results are obtained from a numerical calculation including all mass effects. 
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Model 


■p i dI \ oc 1 + A cosB 
T d cos 6 


MSSM with see-saw mechanism 


A = 1/6 


MSSM with 7?— parity violation: 

"R" (A i2 2A*23 > Aj32A* 22 ) 
"L" (Ai22A* 23 < Aj32A* 22 ) 


A = 1/6 
A = -1/6 


Littlest Higgs model with T~ parity: 

"Y" (KZ « Cdd) 


A = 5/18 
A= 1/2 


Xop colour —assisted technicolour: 

"Z" K»m^) 
"P" {m^m z ,) 


A = -5/14 
A = 


Models with doubly-charged Higgs bosons: 
Higgs triplet model 
Zee-Babu model 


A = 1/2 
A = -1/2 



Table 1: The slope of 1/r dr/dcos0 for the various new physics models discussed in Section El The 
numbers are obtained from the analytic results for massless muons. 

new physics models which predict the decay r — * [i[i]x at rates that may be observable at the LHC: the 
i?-parity-conserving MSSM including see-saw neutrino masses at large tan/3, the -R-parity-violating 
MSSM, the Littlest Higgs model with T-parity, the topcolour-assisted technicolour model, and models 
with doubly-charged Higgs bosons. For these models, our Monte Carlo studies of the decay r — > fx/ip, 
at the LHC predict experimental acceptances of approximately 25-30%. We have emphasized that the 
models can be distinguished from each other by measuring the angle between the r polarization vector 
and the momentum of the anti-muon. This can be achieved within a year of full luminosity at the LHC, 
if the branching ratio for the decay r — > fifip, is close to its current upper bound. 
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